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hyperbolic internal dynamics is presented. The present paper integrates stable inversion techniques (that 
achieve exact-tracking) with approximation techniques (that modify the internal dynamics) to circumvent the 
nonhyperbolicity of the internal dynamics - this nonhyperbolicity is an obstruction to applying presently 
available stable inversion techniques. The theory is developed for nonlinear systems and the method is 
applied to a two-cart with inverted-pendulum example. 
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1 Introduction 


Precision output tracking controllers are needed to meet increasingly stringent performance requirements 
in applications like flexible structures, aircraft guidance, robotics, and manufacturing systems. While ex- 
act tracking of minimum phase systems is relatively easy to achieve through approaches like input-output 
linearization [1], output tracking of nonminimum phase systems tends to be more challenging due to fun- 
damental performance limitations on transient tracking performance [2]. This poor transient performance 
has been mitigated by using pre-actuation in the stable-inversion based approaches [3, 4]. However, the 
preactuation time (during which most of the preactuation control effort is required) depends on the unstable 
poles of the linearized internal dynamics - the preactuation time increases as the unstable poles approach the 
imaginary axis. In the limiting case, with the poles on the imaginary axis (nonhyperbolic internal dynamics), 
presently available inversion-based techniques are not applicable and the effective preactuation time tends 
to become infinite. The present work extends previous results for linear systems in [5] to output tracking of 
nonlinear nonminimum phase systems, which have nonhyperbolic internal dynamics. 

Output tracking has a long history marked by the development of regulator theory for linear systems 
by Francis and Wonham [6] and the generalization to the nonlinear case by Byrnes and Isidori [7]. These 
approaches asymptotically track an output from a class of exosystem-generated outputs. Further, extensions 
to the Byrnes-Isidori regulator have been described in [8]. The main problem with the application of these 
techniques to the output tracking of nonlinear systems is computational. While the linear regulator is 
designed by solving a manageable set of linear equations, the nonlinear regulator design requires the nontrivial 
solution of a first order partial differential algebraic equation [9, 10] . In contrast, inversion-based approaches 
avoid this computational difficulty and trade it to solve the exact tracking problem for a single desired 
output trajectory rather than solve the asymptotic tracking for a class of outputs. Another problem with 
the regulator approach is that the exosystem states are often switched to describe the desired output - this 
leads to transient tracking-errors after the switching instants. Such switching caused transient errors can 
be avoided by using inversion-based approaches to output tracking [3, 11]. Thus, it is advantageous to use 
inversion-based output tracking when precision tracking of a particular output trajectory is required. 

Inversion, which is key to our approach, was restricted to causal inverses of minimum phase systems in 
the early works by Silverman and by Hirschorn (e g., [12, 13]) because these approaches lead to unbounded 
inverses in the nonminimum phase case. Di Benedetto and Lucibello [14] considered the inversion of time 
varying nonminimum phase systems with a choice of the system’s initial conditions. Rather than choose 
initial conditions, preactuation was used by noncausal stable inversion techniques developed in [3, 4, 15]. 
Such noncausal inverses, which require preactuation, have been successfully applied to the output tracking 
of flexible structures [16, 17], and aircraft and air traffic control [4, 18]. There is, however, a fundamental 
limitation to the presently available inversion techniques - they are only applicable if the internal dynamics is 
hyperbolic, and inversion-based output tracking has been a challenge for systems with nonhyperbolic internal 
dynamics. 

Although Huang [8] has proposed some sufficient conditions for developing a regulator for systems with 
nonhyperbolic internal dynamics, the regulator design remains computationally difficult. There are also 
several approximation based output tracking techniques, where the central philosophy is to replace the 
internal dynamics with a dynamics that provides satisfactory behavior, and then to develop the controller 
based on the altered system [19, 20, 21]. The techniques most relevant to this paper are developed by 
Gurumoorthy and Sanders [19], and by Gopalswamy and Hedrick [21] - output redefinition (or modification) 
is used in these to approximate the unstable internal dynamics with a stable system. Such approximation 
based approaches are integrated, in this paper, with stable-inversion based techniques to achieve inversion 
of linear systems with nonhyperbolic internal dynamics in [5]. However, the present technique does not use 
output modification to stabilize the internal dynamics, rather output modification is only used to remove 
the nonhyperbolicity. Additionally, for near nonhyperbolic systems, the present approach allows a tradeoff 
between the precision-tracking and the amount, of pre- actuation- time that is needed in the control effort. 
This tradeoff between stable-inversion and approximation of the unacceptable internal-dynamics has been 
studied for helicopter-hover control in [5]. The present work extends the results for linear-systems in [5] to 
nonlinear systems. 

We begin with a background on the stable-inversion based output tracking technique, and computational 
issues are presented in Section 3. The technique is applied to an example in Section 4, and conclusions are 
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in Section 5. 


2 The Stable Inversion problem 


2.1 Inversion-Based Output Tracking Scheme 

Here we describe how the inversion approach is used to develop output tracking controllers. Consider a 
system described by 

z(t) = f[x(t)] +< 7 [x(t)] u(t) 

( 1 ) 

y(t) - A [*(<)] 

where y(f) = [j/i (<), xj 2 {t ), . . . , y P (f)] T the output, with the same number of inputs and outputs, i.e., 
^ ^ ) an< ^ x(t) G 3? n is the state. The functions /(*),y(*) and /i(-) are assumed to be sufficiently 
smooth with /( 0) = 0 (i.e., x — 0 is an equilibrium point) and A(0) = 0. 

Let yd{-) be the desired output trajectory to be tracked. In the inversion-based approach we, first, find 
a bounded, nominal input-state trajectory, [uyy(-), x re f(-)) that satisfies the system equations (1) and yields 
the desired output exactly, i.e., 


Xrej{t) — f[ x ref(t)] + 9 [ x rcf ( 0 ] u Jf (0 
Vd{t) ~ h[x re j(t)} 


( 2 ) 


and, second, we use the exact-output yielding input trajectory, ti^(-), as feedforward and the system is 
stabilized by using feedback. It is noted that in this output tracking scheme, the feedforward input ujj ( •) is 
computed off-line. 

In the absence of modeling errors, initial condition errors, and external disturbances, exact-tracking can 
be achieved by the use of feedforward alone. If the system (1) is stable (or stabilized with feedback before 
the inversion is applied) then standard Lyapunov arguments [22] can be used to show that, with the use of 
feedforward, x{t) -4 x rcf {t) and y{t) -4 y d (t) as t -4 00 and output tracking is achieved [3, 15]. The feedback 
scheme can also include integrals of the tracking error [22] to obtain zero steady-state errors (see Figure 1). 
While stabilization of the system can be achieved through standard techniques, the main challenge is to find 
the inverse input-state trajectory [ujj (•), x r<! /(-)] - especially for systems with nonminimum phase dynamics. 
This paper addresses the issue of generating the inverse input-state trajectory for a system and can be used 
in conjunction with any feedback stabilization scheme. 



Figure 1. The Control Scheme 


2.2 The Internal Dynamics 

In this subsection, it is shown that finding the inverse input-state trajectory is equivalent to finding bounded 
solutions to the system’s internal dynamics. 


Assumption 1 System (l) has a well defined vector relative degree , r [rj, r 2 , . . r p ]. 

The well-defined vector relative degree assumption enables the system equations to be rewritten, through a 
co-ordinate transformation 


*(0 = f(x(t)) 


' C(t) 

. v(t) 


ncwM*)) 


m ■ 
m J 


( 3 ) 
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in the following form [1, 21] 



Cl(<) = ^CiCl(<) + ^CaC2 (<) 



6(0 = si[C(i).v(0] + *a[C(0> , K0M < ) 

(0 


n(t) = *3[C(0.»;(0] + 


where £ represents the output along with its time derivatives, £ is partitioned as | 
given by 

1 T 

Ci(<) i Ca(t) J and 15 

±L i «J 3 !/i 

L yi at Ho ■ ■ 

<T*" 3 yi _. d r *~ 2 .. d r p “ 3 .. : d r I -1 

* dt r I - 3 y* *•’ ■■■ yi dt^-^y^ 

< 5 > 


Further A ^ is a block-diagonal matrix, A ^ = diag[Aj , A 3 , . . .,A P ], with each block consisting of zeros 

except for ones on the super-diagonal (A* is a (r* - 1) x (r* — 1 ) matrix for each 1 < k < p). A ^ is a matrix 
whose elements are zeros except for elements on the (rj +r 2 + . . . -f r* — k)-th row and fc-th column, which 
are ones (for all 1 < k < p). Further, from Assumption (1), •) is invertible in a neighborhood of the 

origin, and since the origin is assumed to be an equilibrium point, we also have si( 0 , 0) = 0,53(0,0) =0. 
Note that the desired £(•) is known when the desired output trajectory t/<*(*) and its time derivatives are 
specified. This desired £(•) is defined as £*(♦). If exact output tracking is achieved (i.e., if £(f) = Q(t)) then 
the control law for maintaining exact tracking can be written as, from equation ( 4 ), 

«(0 = MWO.'HOr 1 [CmW - s,[Cd(<).v(<)]] (6) 

which results in state-equations of the form 

c(o = 6 (<) ( 7 ) 

v(t) = s 3 [c(t) } Tj(t)] + [<(*), >?(0M0 

= «»[c(0, «?(<)] + *i(0] muo, m)r l [cw(o-*i[c-(o. *?(<)]] 

( 8 ) 

This is the inverse system and equation (8) is referred to as the internal dynamics. Solving the internal 
dynamics is key to finding the inverse input state trajectories. If a bounded solution, 77 <*(■), to the internal 
dynamics (8) can be found, then the feedforward input can be found through equation (6) as 

«//(<) = (*2(c4<).^(<))r' [cj.d(o-*i(c-(<), %(<))] ( 9 ) 

and the reference state trajectory can be found as 


x rcf(t) — jF[C<f(0i )] 


C <(0 


( 10 ) 


Thus, finding a bounded solution, Tjd(-), to the internal dynamics (Equation 8) is key to finding the inverse 
(bounded) input-state trajectories (Equations 9 and 10) which are needed to implement the inversion-based 
output-tracking scheme shown in Figure 1. 


2.3 Modified Internal Dynamics 

Standard inversion schemes [12, 13] that integrate (forward in time) the internal dynamics (8) lead to un- 
bounded solutions if the origin of the internal dynamics is unstable (nonminimum phase systems). Noncausal 
inversion (e.g., [3]) leads to bounded but noncausal solution of the internal dynamics. While significant im- 
provement in output tracking performance is possible, such stable-inversion techniques are not applicable to 


3 


systems with nonhyperbolic internal dynamics. In this subsection a compromise between stable inversion and 
approximation-based inversion schemes is proposed. The key is to modify the internal dynamics by giving 
up exact output tracking - enough to remove the nonhyperbolicity, and then to apply stable-inversion. The 
difference between the proposed technique and other approximate techniques is that the internal dynamics 
is perturbed only to remove the nonhyperbolicity, and not to stabilize the internal dynamics as in other 
approximation schemes [19, 21]. 

To modify the inverse system an extra term, is added to the control law (6) as follows 


u (<) = [ S 2 [C(<), »?(<)]] 1 [C 2 ,d (0 -*i[C( 0 . *?(<)] + «(*)] 


( 11 ) 


Note that the modified input does not attempt to achieve exact tracking of the desired output (i.e., £(*) = Q(.) 
is not required) Therefore, the modification, u(<), of the control law can be used to modify the internal 
dynamics. With this control law, the system equation (4) becomes 

Ci(<) = AMt) + A^t) 


6(0 = 6 , d (0 + v (0 


n(i) = s 3 m,v(t)} + s 4 [at),v(t))[s 2 [C(t),v(t)]]~ l [ 6 , 40 -si [C( 0 ,*K 0 ] + *(<)] 

For ease in notation, we define the tracking error, as 


e ( (<) := 


’ e c>(<) ' 


. e Ca(0 . 

. 


C>(<)-Cm(*) 1 

C 2 (<) - C 2 ,d(t) J 


( 12 ) 


and rewrite Equation (12) as 


«C. (<) 

= ^Ci e Ci(*) + 


«(»(*) 

= »(0 


W) 

= s 3[e ( (0 +Ci(<), »?(<)) + 

(13) 


s 4 [e((<) + Cd(<). * 7 ( 0 ] [ s 2 [ec (0 + 0 (<), »?(^)]] ' [C 2 ,d (0 - «i [«<(<) + Q(<)> »?(*)] + »(<)] 


: = s( e c(*).»?(0> y 4<). *>(<)] 


where 


Y d (t) 


■ 0(0 ' 
. 6,<f(0 . 


If t>(<) = 0 then Equation (13) represents the inverse-system, which has to be solved to find the exact- 
tracking inverse input-state trajectories. However, the internal dynamics part of equation (13), 


n(t) = i[e c (0,»?(0.y<<(0,o] , 

is nonhyperbolic and may not have a bounded solution (see the Example in Section 4) - thus the inversion 
problem may not have a bounded solution. Even when bounded solutions exist for the nonhyperbolic internal 
dynamics, currently available computational techniques for stable-inversion cannot be directly used since they 
require hyperbolicity of the internal dynamics [3, 4]. Our approach is to appropriately choose v( ) to remove 
the nonhyperbolicity of the inverse system (13), which includes the internal dynamics. We begin with the 
following assumption. 


Assumption 2 System (1) is controllable tn the first approximation [1], 

The above assumption implies that the modified inverse system is also controllable in the first approximation 
since the difference between the original system (1) and the modified internal dynamics (13) is only a 



co-ordinate transformation and a static state- feed back (see [21] for a similar argument used for output 
redefinition). Next, a feedback of the form 

•o = f [ ] < i4) 

is chosen such that the modified inverse system is hyperbolic - i.e., all poles on the imaginary axis are 
moved. Note that this change to a hyperbolic system can be achieved through an arbitrarily small F since 
nonhyperbolicity is not a structurally stable property. With this control law Equation (13) becomes 


£C. (0 = 4(,e (l (<) + A ( ,e (j (t) 
«<*(*) = F 


e ( (t) 

l(0 


i(0 = • J 

which is re-written in a simplified form as 

T(t) ] = 

Next, stable inversion of the modified inverse system (16) is carried out [3]* 


(15) 


(16) 


3 Computation of the Inverse 

This section discusses the computation of the inverse for the hyperbolic, modified, inverse-system. First, an 
iterative algorithm is presented - one of the steps involves finding bounded solutions to a linear (unstable) 
system. Second, the explicit technique to find bounded solutions is presented, and the amount of preactuation 
required for implementing the inversion-based controller is discussed. We begin with the algorithm to find 
bounded solutions to the modified internal dynamics [3]. 


3.1 Iterative Algorithm to find Inverse 


• Step 1 

Rewrite the modified inverse (16) as 


± r e f (<) 

dt [ r)[t) 



+ 


|s(e f (i), »?(<), y d (<)) 



(17) 


where the first term on the right hand side (r.h.s.) is the linearization of s(*,-, ■) with respect to the 
first two variables, e^( ) and tj(-), and the second term on the r.h.s. represents a perturbation. This 
motivates the following iterative solution of the internal dynamics. 


• Step 2 

We begin with [e^(*) ^(<)]^ = 0 and then, at at each iteration in the following scheme, the bounded 
solution of the following linearized system-equation is found. 


d 

Jt 


ec(*) 

»?(<) 


J jV+1 


e C (0 ' 

v(t) . 


+ 


7V + 1 


js ( e (f l,N){t),VN(t),Y d (t)) 


e<(<) 1 \ 

l(0 J N l 


where 


' e<(<) ' 


e c,/v+i(0 

. 1(0 . 

N+l 

VN + l(t) 


and N = 1, 2, 3 


(18) 

( 19 ) 
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3.2 Bounded Solution to Unstable Linear System 

In the above iterations, solving equation (18) implies finding a bounded solution to a linear, hyperbolic (but 
potentially unstable) system - this is described next. 

We begin by decoupling the linear system (18) into stable ( z ,) and unstable (z u ) subsystems. Since the 
modified internal dynamics is hyperbolic, there exits a decoupling transformation U such that equation (18) 
can be re-written as 

z (u,N+l)(0 = Su*(ti f AT+l)(f) +G\ i Yw(<) 


where 

*N+.(<) := [ ! (, ' JV+,) S 1 := U \ ] . 

L Z (u,N+l)(<) J [ T) N+i (t) J 

and 

An approach to find bounded solutions is to enforce the boundary conditions that Z( s ,n+ i)(-oo) = 0 and 
z (u,tv+ i)(°°) = 0. This leads to unique bounded solutions by flowing the stable subsystem forward in time 
and flowing the unstable system backward in time - this yields (see [15] for computational issues) 

*(.,W+l)W = floc eS - (t - T)G ’ Y "(T)dT Vt€ (-00,00), 

, (20) 

z (u ,N+i)(t) = £ 0 e- s ’‘WG u Y N (T)dT V/e (-oo.oo). 

Next, a change of co-ordinates yields the bounded solution to (18) as 

c <(*) 1 -[/-if *(»,n+i)( 0 l . . 

9(0 Jn + 1 l^ + i)(Oj N+1 - (21) 


*N+.(<) := l ( 
Z ( 


3.3 Convergence of Iterative Algorithm 

The following Theorem states that, as N — > oo, the solutions of the above iterative algorithm converge to a 
bounded solution of the modified inverse system (16). 


Theorem Let the hyperbolic, modified, internal dynamics s(-, •, •) and its linearization S satisfy the fol- 
lowing Lipschitz like condition on a neighborhood of the origin, (0,0,0), 


*(«i,»7i. Yi) - S 


]}-{ 


s { e 2, Y 2 ) “ S 


+ K 2 II Yi - Y 2 


Then, provided A 1 , A' 2 and || Y<*(*) are small, the iteration scheme converges to a solution of the modified 
inverse system (16), i.e, 

e C e «,<0 1 

. V'In % J 


as N — > oc, and satisfies the modified 

Vd 


inverse system. 


Proof 

See [3, 4]. 


Thus, stable inversion technique is applied after removing the nonhyperbolicity in the internal dynamics. 
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Remark 1 The conditions in the theorem are sufficient conditions for finding bounded solutions to the 
internal dynamics . They are not , however, necessary for the existence of bounded solutions to the internal 
dynamics - bounded solutions can exist even when the internal dynamics are non-hyperbolic . Further , S in 
the above theorem can also be chosen to be different from the linearization of the nonlinearity s, however, 
additional conditions have to be imposed (see, [3], Condition 1 in Definition 5). 

This completes the inversion-technique for nonminimum phase systems with nonhyperbolic internal dy- 
namics. To summarize, the reference state trajectory is found as (from Equation 10) 


x ref{t) — , 7)<i(t)] 


’ C-(0 
. w(t) 


where 

0(0 :=: 0(0 + e ( c ,<*)(0 

and the feedforward input, «//(•), is found from equation (11) and equation (14) as 

«//(*) = [*2[Cd(0, w(<)]] -1 [c 2 ,*(<) - «i[Cd(0, •?(*)] + F [ ]] (23) 

The inverse input-state trajecotry, [ti //(•)> z re/(')]i i s then used with feedback stabilization to obtain output 
tracking [3j. 


3.4 Error Analysis 

The following Lemma establishes an upper bound on the tracking error, for a given choice of feedback, 

F, that removes the nonhyperbolicity of the inverse system (15). 

Lemma 1 Let 


= U~ ] 


I(t)e s,t 0 

0 — /(— 0 e ” 5uf 


where /(•) is the unit step function, 5 is defined from equation (17) as 


5 = 



/t Ca 0 
F 
S 


where 5* is the linearization of s(-, •, •) with respect to the variables, e^(-) and f]( ) (see equation (15)), and 
S 9 and S u represent the decomposition of linearized inverse system into stable and unstable parts by the 
co-ordinate transformation, U, as discussed in Subsection 3.2. Further, let 


II Mill = nmax[|^(.)|| r 

Then the tracking error for a particular choice of F can be bounded as 


ll*C<Wlleo < 




Proof This follows from Theorem 3.6 in [23]. 


3.5 Preactuation Time 

Stable-inversion techniques overcome fundamental limitations on transient tracking performance of non- 
minimum phase systems [2] by using preactuation. However, the preactuation time (i.e. } when most of the 
preactuation effort is required) tends to be unacceptably large if the unstable poles of the linearized, modified 
internal dynamics are close to the imaginary axis - this dependence is established in the following Lemma. 


i I 
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Lemma 2 
Let 


• all the unstable poles of linearized, modified inverse system (eigenvalues of S u ) lie to the right, in the 
complex plane, of the line Real(s) = a for some positive a, and 


• the support of Yd(-) lie in [J 0 , oo) for some to. 


Then there exists a positive scalar, M, such that the bounded solution to the internal dynamics (defined by 
the Theorem) satisfies 


e ((,< t )(0 

» 740 




for all time, t , before the start of the maneuver at <o- 


Proof Since the maneuver starts at /o 5 any solution to the modified internal dynamics (16) satisfies the 
autonomous equation, 

= s ( e c(0. *?(<). 0)- (24) 

for all t < t 0 . Further, any bounded solution to the autonomous equation must lie on the unstable manifold 
before the start of the maneuver. The rate of convergence to zero as time tends to -oo, of solutions that 
lie on the unstable manifold of the nonlinear modified inverse system, directly depends on the location of 
poles of the linearized dynamics. In particular, the rate of convergence depends on the unstable poles of 
the linearization (which are the eigenvalues of S u ). The existence of a positive scalar M that satisfies the 
statement of the lemma follows from the saddle-point property (see [24] Theorem 6.1, and [251 Theorem 
19.9). a 

The present technique also provides a way to reduce the preactuation time for systems with near nonhy- 
perbolic internal dynamics. The Lemma states that the desired state-trajectory exponentially tends to zero 
as we go back in time before the start of the maneuver at t 0 . Then the continuity of the input with respect 
to the state (smoothness of /, g and h, and the well defined relative degree assumption) implies that the 
preactuation input also tends to zero exponentially. The rate at which the preactuation becomes zero can 
be increased by moving the unstable poles of the linearized, modified inverse system (i.e., the eigenvalues 
of S u ) away from the imaginary axis - by appropriately choosing F in equation (14). This reduction in 
preactuation-effort is obtained at the expense of exact output tracking (see [5] for a linear example). 
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Figure 2. Example: Two Carts and Inverted Pendulum 


Here, the inversion-technique is applied to an example two-cart and pendulum system shown in Figure 
2. The inverted pendulum on a cart has been well studied in literature (see, for example, [19]), and has 
nonminimum phase dynamics - the internal dynamics is hyperbolic. An extra cart, is added here, which 
introduces nonhyperbolicity in the internal dynamics. The input, to the system is the applied force, F(t), 



and the output is the position, x\ (f), of the cart carrying the inverted pendulum (see Figure 2 ). The equation 
of motion for the system can be obtained as 


(M + m)xi (t) + ml cos 0(f) 0(f) 
m/ cos 0(f) xj(f)-{-m/ 2 0(f) 
Mx 2 {t) 

which can be rewritten in state-space form as 


F(t) + mld 2 (t) sin 0 (t) — K (xi(<) — x 2 (t)) 
mgl sin 6 (i) 

-K (x 2 {t) - 1,(0) 


(25) 


Cl (0 = 

6(0 

[ m/m 



i 

m = 

■»(*>») 

+ lrj 4 (t ) 2 sin r) 3 (t) - K/m(£i(t) - Tji(t)) 

- <7Cos^ 3 (<)sin 77 3 (<)j 

m{t) = 

mit) 



V2 (<) = 

(K/M) (Ci(<) • 

-»?i(0) 



m (<) = 

m{t) 





»?4(0 = 

1 I 
Moiy | 

[~»/4(0 2 

sin 773(0 cos 

%(0 



^ mT 

{{M + m)g sin t? 3 (7) 

- F(t) cos 7 ^ 3 (<) + K (Cl (0 - 

7,i(<))cOS T) 3 (t)}] 


where Cl :=*i, C2 *1, m := x 2i 7 ft := x 2l rj 3 := 0, and r } 4 0, and 7(773(0) := M/m + [sin 773(f)] 2 

Given the desired output trajectory profile, the input that maintains exact tracking (i.e., which maintains 
C (0 — Cd( 0 ) obtained from equation ( 26 ) as 

F(t) = my{r) 3 )x lid (t) - m/[r; 4 (f)] 2 sin rf 3 (t) + K [Ci (f) - (f)] + mg[cos 773(f)] sin 773(0 


With this exact tracking control law, the linearized internal dynamics (represented by 77), is given by 



m (0 ’ 


0 

1 

0 

0 ' 


»M0 " 


0 

0 ' 


d 

V2 (0 


-K/M 

0 

0 

0 


*? 2(0 

4 - 

K/M 

0 

r *i,d(o 

dt 

173(7) 


0 

0 

0 

1 


173(0 

0 

0 

L *m (0 _ 


. »? 4(0 . 


0 

0 

9/1 

0 


. 174(0 _ 


0 

-l/l _ 



Note that the linearization has two poles at ±y/—K/M on the imaginary axis (in the complex plane) which 
correspond to the cart dynamics (x 2s x 2 ) with xj fixed; i.e., a spring-mass system. The other two poles are 
at ±^/gTl which correspond to the linearization of the inverted pendulum dynamics ( 0 , 0 ). The two poles 
on the imaginary axis lead to nonhyperboiic behavior and stable-inversion techniques fail - this corresponds 
to the case which requires infinite pre-actuation. 


Remark 2 Note that the nonhyperboiic spring-mass system will not have a bounded solution if the desired 
output excites the spring-mass system at its resonance frequency . Thus the nonhyperboheity of the internal 
dynamics can imply that exact-tracking cannot be achieved for some desired- output trajectories. 

By modifying the desired trajectory, the internal dynamics can be approximated by a hyperbolic system, 
for which existing stable-inversion techniques are applicable. This modification is achieved by adding an 
extra term ^(*) to the control as follows 


u //(0 - ^(0 = ™ 70 fe)[iM (0 + v (0] — r/z/ [774 (^)] 2 sin 

+# [Cl io - * 7 i( 0 ] + "i0[cos 773(f)] sin 773(f) 

The modified inverse system becomes 


^c. (0 

= C (,(0 

0 

= *(0 

* 7 i (0 

= 172(7) 

1)2(0 

= —K/Mrjx + I</M[x 1 >d (t) + e Cl ( 0 ] 

1)3(7) 

= i 74 (<) 

1)4(0 

= 9/1 sin 7/3(7) - 1//cost? 3 (0[£m( 7) + v(7)] 


( 29 ) 



where 


V' 


v 


v 


V 






e ( (t) := 


At) 1 

Ci(t)-x M (<) 1 

tUi) J 

. C 2 (f) - i M (0 J 




Figure 3* Results of Inversion 

The additional input, v, is chosen (v = F[eJ r] T ] T ) to remove the nonhyperbolicity. This can be achieved, 
for example, by pole-placement algorithms. Simulation results are presented next. The system parameters 
were chosen, for the simulation, as A' = 10 N/m, M = 1 Kg, l = 9.8m, g = 9.8 m/s 2 and the F used in 
simulations is 

F = [— 4.52e — 2 - 3.04e - 1 5e - 4 1.5e-3 0 0]. 

This F removes the nonhyperbolicity of the internal dynamics and stable inversion of the modified inverse 
system is carried out using the algorithm in Section 3. Simulation results are presented in Figure 3, which 
shows the error caused in output-tracking due to redefinition of the output trajectory. The tracking error is 
less than 1% of the maximum value of the desired output (see Figure 3) - the rate at which the error goes 
to zero depends on the choice of F. The internal dynamics ( x 2 ,0 ) is also shown in Figure 3. Note that 
the nonhyperbolicity is circumvented and stable inversion is achieved with a relatively minor modification 
of the output . 

5 Conclusion 

A technique to achieve output tracking for nonminimum phase nonlinear systems with nonhyperbolic internal 
dynamics was presented. The approach integrated stable inversion techniques (that achieve exact tracking) 
with approximation approaches (that modify the internal dynamics) to remove the nonhyperbolicity of the 
internal dynamics. It was shown that, by giving up some of the precision in tracking, it is possible to achieve 
stable inversion of nonlinear nonminimum phase systems with nonhyperbolic internal dynamics. 
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